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Abstract. In this paper we emphasize the close connection between the chemical and spectrophotometric evolution 
of stellar systems: Chemical yields from stars correspond to a precise fraction of their emitted light. We translate 
this concept quantitatively. Starting from simple stellar populations, we derive useful analytical relations to 
calculate the stellar fuel consumption (emitted light) as a function of basic quantities predicted by stellar models, 
i.e. the mass of the core and the chemical composition of the envelope. The final formulas explicate the relation 
between integrated light contribution (total or limited to particular evolutionary phases), chemical yields and 
stellar remnants. We test their accuracy in the case of low- and intermediate-mass stars, and indicate the way to 
extend the analysis to massive stars. This formalism provides an easy tool to check the internal consistency between 
the different stellar inputs adopted in galaxy models: The fuel computed by means of the analytical formulas 
(corresponding to a given set of chemical yields) should be compared to the exact values given by the luminosity 
integration along the stellar evolutionary tracks or isochrones (corresponding to a given set of spectrophotometric 
models). Only if both estimates of the fuel are similar, the stellar inputs can be considered self-consistent in 
terms of their energetics. This sets an important requirement to galaxy models, also in consideration of the fact 
that different sources of input stellar data are frequently used to model their spectro-photometric and chemical 
evolution. 
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1. Introduction chemical yields for CEMs. Other distinctive featm^es of 

Ci i,n these fields are that in SPMs star formation rates and 

Galaxy evolution models represent fundamental tools for ^ n- -^ i i- ^^ ■ ^ ^ ^^ i i 

,. . _ ... „, -11 age-metallicity relations are usually inputs to the models, 

the interpretation of a wide variety of observational data, , . ^t^tv r J^ , , r .^ i i ; ■ 

whereas m OEMs they are outputs of the calculations, 
going from star counts m our Galaxy to galaxy counts 

in deep fields. One of the essential assumptions in these This picture is nowadays rapidly changing, since the 

models is that we know how stars convert their nuclear present trend is to model complex stellar aggregates like 

fuel into light and newly- synthesized elements. Then, this galaxies accounting simultaneously for all their aspects, 

knowledge is used to model the light emission from galax- mainly spectrophotometric and chemical, and possibly 

ies, as well as the chemical enrichment of the galaxy including the dynamical one. Some examples of unified 

medium and the baryonic mass locked into stellar rem- SPMs and CEMs of galaxies can be found in Arimoto & 

nants. The basic formalism is fully described in the pio- Yoshii (1986, 1987), Bressan et al. (1994), Einsel et al. 

neering work by Tinsley (1980). (1995), Vazdekis et al. (1996), Chiosi et al. (1998), and 

However, it is also a fact that light and chemical evolu- Ferreras & Silk (2000). The simplest way to couple GEM 

tion of galaxies have been historically developed into two and SPM is quite obvious: one may simply re-direct the 

separate disciplines, that we refer to by the generic names output from GEM as input to SPM and, as a matter of 

of "spectro-photometric models" (SPM) and "chemical fact, this choice has been frequently made in literature. In 

evolution models" (GEM) of galaxies. This separation has this way the models naturally succeed in reproducing some 

probably been the result of the very different basic tools observed scaling relations, like for instance the colour- 

and input data required in these two fields, namely: stellar magnitude (mass-metallicity) relation of ellipticals (e.g. 

theoretical isochrones and spectra for SPMs, and stellar Bressan et al 1994). 
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cent indications by Ibata et al. (2000; see also Chabrier 
1999) that a significant fraction of the detected massive 
compact halo objects (MACHO) could consist of white 
dwarfs, would pose non-trivial problems to present chcmo- 
spectrophotometric models of galaxies: high rates of white 
dwarf formation at early epochs would not be compati- 
ble with the observed deep galaxy counts (Chariot & Silk 
1996), and with the modest chemical enrichment now ob- 
served (Fields et al. 2000). Clearly, these studies point at 
a link between the light of distant galaxies, and chemi- 
cal abundances and stellar remnants in the local universe. 
In this regard, we should also mention the development 
of more "global" analyses (e.g. Madau & Pozzetti 2000; 
Fall 2001), that consider the constraints provided by the 
cumulative emission of cosmic structures (the optical to 
far- infrared extragalactic background light), and the av- 
erage metallicity and stellar mass density of the present 
universe. 

However, the fact that nowadays different aspects of 
galaxy evolution are being considered simultaneously, does 
not necessarily imply that the resulting models are self- 
consistent as a whole. A potential source of inconsistency 
could reside in the (quite frequent) use of different sources 
of stellar data as input to SPM (e.g. luminosities and 
lifetimes) and CEM (e.g. chemical yields and remnant 
masses). 

In this context, the present paper emphasizes the close 
relation between the chemical and spectrophotometric 
evolution of galaxies. We derive the expected relations 
between the integrated light of a stellar population, the 
chemical yields of their component stars, and the mass 
of stellar remnants. These provide important constraints 
to any model of galaxy evolution, which have so far been 
either not explicitly checked for, or simply neglected. 

Wc propose such a consistency check may be carried 
out with the aid of an analytical formalism, based on sim- 
ple relations involving stellar parameters. The check can 
be performed as follows. First one derives the fuel (emitted 
light) directly from the luminosities of the adopted stellar 
tracks or isochrones (Sect, g), and then compares it with 
the fuel analytically derived (Sect, y), which corresponds 
to the particular set of chemical yields in use. Of course, 
the test is successful if the two determinations converge to 
the same result, as fully discussed in Sect. ^. 

2. Single-burst stellar populations 

Single-burst stellar populations (SSPs) are the building 
blocks of the composite populations we call galaxies. The 
integrated stellar light in the pass-band A of an SSP of age 
t is given by Q 



(1) 



Lf''{t)= / cj){Mi)Lxt{Mi)dMi 

JQ 

where M; corresponds to the initial stellar mass, L\t{Mi) 
denotes the luminosity along the isochrone of age t, 

^ Effects such as extinction, gas emission, etc. can be consid- 
ered a posteriori in the models. 



and 0(Mi) is the initial mass function (IMF). This lat- 
ter is usually normalized so that the total SSP mass is 
equal to a known quantity Mt, e.g. Jp M; (l>{Mi) dMj — 
Mt. Integrated magnitudes and colours follow straightfor- 
wardly from the quantities L^^{t). 

From the above equation, and taking a number of jus- 
tified approximations (see Girardi & Bertelli 1998), we can 
also recover the so-called "/jtel consumption theorem!'' of 
Renzini & Buzzoni (1986): 



Kint) - Km + ^H h{t) > F,(AfT0) , 



(2) 



In brief, this equation tells us that the integrated bolo- 
metric luminosity of a SSP of given age can be expressed 
as the sum of the integrated luminosity along the main se- 
quence (MS), ibcl(Oi ^'^d ^ quantity that is proportional 
to the total fuel, Ft = ^ ■ Fj{Mt^o), consumed during all 
the post-main sequence phases of a star with mass equal 
to the turn-off mass, Mto ■ In this equation b{t) is the evo- 
lutionary flux, i.e. the number of stars leaving the MS per 
unit time. The fuel consumption, Fj, represents the bolo- 
metric light emitted by all stars in the post-MS evolution- 
ary stage j. However, instead of being expressed in energy 
units, this light is expressed in units of "equivalent-mass 
of hydrogen burnt" (e.g. in Mq). This energy-to-mass con- 
version is properly taken into account by the constant Ah 
in Eq. (g), which expresses the efficiency (i.e. energy re- 
leased per unit mass) of H-burning reactions. In this work 
we adopt Ah = 9.75 x 10^" L© yrM©"^, which is derived 
from the net Q-value of H-burning reactions via the CNO 
cycle. 

Since the burning of a given mass of helium provides 
only one tenth of the energy that is provided by the 
nuclear consumption of the same mass of hydrogen (i.e. 
Aho ^ 0.1 Ah), Fj{Mto) is then given by 

Fj {Mto) - AMh^ + 0.1 AAf hc, , (3) 

where AMhj and AA/hoj are the masses of H and He, 
respectively, nuclearly burnt during the j-th post-main se- 
quence evolutionary phase of the star of initial mass Mto- 
From Eqs. (w) and (0) above, it is also clear that 



F,(Mto) = 



1 



Ah b{t) 



^{Mi)Lt{Mi)dMi, 



(4) 



where the integration is limited to the section of the 
isochrone with age t that corresponds to the j'^- 
evolutionary stage only. 

On the other hand, if we are dealing with stellar evolu- 
tionary tracks instead of isochrones, the fuel of a star with 
initial mass Mi can be calculated as the time-integral of 
the luminosity over the duration of that phase, divided by 
the constant Ah: 

F,{Mi) = ^lLM,{t)dt. (5) 

In the previous equations, Fj is expressed in Mq, L is 
given in Lq, and the independent variable t is expressed 
in yr. 
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From the aforementioned definitions it is immediately 
clear that the Fj are important quantities, because they 
relate basic properties of stellar structures (Eq. 0) to the 
integrated light of SSPs (Eq. Eh. Examples of the practical 
application of the above concepts can be found in e.g. 
Greggio & Renzini (1990) and Renzini (1998). It is just 
the stellar fuel consumption and its main features the main 
subject of the analysis we will develop in the following. 



3. Computing the fuel consumption 

Let us now concentrate on the problem of determining, in 
a generalised fashion, the quantities AMhj and AMhgj, 
or equivalently Fj, corresponding to a stellar model of 
given initial mass M^. 

We can devise two methods respectively based on 

1. "exact" numerical integration along the track 

2. analytical approximations 

With the first approach the fuel consumed during the 
j*'^-evolutionary phase can be calculated from the stellar 
track by simply using Eq. (||), or, alternatively, by inte- 
grating along the corresponding section of the isochrone 
with Mto = Ml as in Eq. (|). 

The second method relies on simple approximate rela- 
tions, involving essentially the size of the fuel-exhausted 
core, and the increase/decrease in the mean mass- 
coordinate of the overlying burning shell. In the context 
of this approach, we will present here an analytical for- 
malism to derive the fuel, developed in detail for low- and 
intermediate-mass stars, and briefly sketched for massive 
stars. The relations are of easy practical application as 
they contain quantities that can be extracted from avail- 
able stellar models. 

If the first approach is to be preferred in terms of accu- 
racy - as it straightforwardly follows from the definition of 
fuel consumption itself (see Sect.0) -, the second approach 
has the merit of expliciting intimate inter-dependences be- 
tween (apparently) different properties of stellar popula- 
tions (i.e. light emission and chemical enrichment). The 
accuracy of the formalism will be tested by comparison to 
the results obtained with the "exact" method. 



3.1. Basic considerations 

Two fundamental ingredients in the calculation of the 
fuel are the mass-coordinates of the H- and He-exhausted 
cores. This is already suggested by Eq. (^), that expresses 
the fuel in terms of two contributions - AA/h and AA/hc 
- which are obviously related to the growth of the H- and 
He-exhausted cores, respectively. In this sense we may say 
that, as the star evolves, the He and CO cores keep track 
of the fuel being consumed, by accreting the ashes of the 
nuclear burnings occurring above them. We may also ad- 
vance that, in general, the derivation of the fuel will in- 
volve the difference between the core masses at some "ini- 
tial" and "final" stages of the phase under consideration. 



From Eq. (H), we can also derive some general rules 
that are valid for any increase of the core mass AM^. If 
we are considering a H-exhausted core that has grown in 
mass by simply transforming H into "^He, the correspond- 
ing fuel is increased by AFj — X AMc, where X is the 
original abundance of H (in mass fraction) in the mate- 
rial nuclearly burnt. If instead, we are dealing with a He- 
exhausted core that accretes the products of He-burning 
(mainly carbon and oxygen via the a-capture reactions), 
we would have AFj =0.1 AA/c- Finally, if during the same 
evolutionary stage, a core has grown in mass through the 
successive conversions of H into He and then He into CO, 
we would have AFj = {X + 0.1) AAfc- 

However, besides these general and quite simple con- 
siderations we should add others which are equally im- 
portant, though (probably) somewhat less intuitive. We 
remark, in fact, that the knowledge of the M^ is not suffi- 
cient to estimate the fuel, since part of the burnt material 
initially deposited in the core may be reduced by dredge- 
up events. The net effect is that products of nuclear burn- 
ings previously occurred in situ in core regions (burning 
shells included) are carried away to the outermost layers 
(dredge- up), and a fraction of them can also be irreversibly 
lost by the star (mass loss). Then, we can already expect 
that this bit of information subtracted from the core is 
intimately related to the stellar yields, as we will better 
explicit in Sect. 3.2. 

To allow an easier understanding of the formalism de- 
veloped in this study let us recall now a few relevant as- 
pects of dredge-up and mass-loss processes characterising 
the evolution of low- and intermediate-mass stars. 

The dredge-up episodes - caused by the penetration 
of the convective envelope into regions previously affected 
by nucleosynthetic processes - typically occur in red-giant 
stages (i.e. RGB and AGB) between major phases of core 
nuclear burning, possibly producing two effects: 

— changes in the chemical composition of the envelope, 
and 

— inward shift of (either one or both) the hydrogen- 
helium and helium-CO discontinuities (or chemical 
profiles) . 

Evolutionary stellar calculations (e.g. Girardi et al. 
2000; Marigo et al. 1999) predict that the I''* dredge-up 
occurs at the base of the RGB in low- and intermediate- 
mass stars of any mass, the 2"^^ dredge-up takes place only 
in stars with initial mass Mi ^3 — 4 Af© at the beginning 
of the Early- AGB (E-AGB), and the 3'''^ dredge-up0 may 
occur at thermal pulses during the TP-AGB phase. 

The 1*** and 2"^^ dredge-up bring up to the surface ma- 
terial which has undergone H-burning, with a consequent 
net enrichment of ^He in the envelope. The third dredge- 
up can contribute to increase the surface "^He abundance 
as well, at the same time injecting into the envelope prod- 



^ this term comprehends all recurrent dredge-up events dur- 
ing the TP-AGB phase 
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Low-mass star 

1. 




WD 



TAMS 



RGB-tip 



Ist-TP 



AGB-tip 



evolutionary stage 

Fig. 1. Schematic representation of the structure of a low-mass stars at several evolutionary stages, namely (from left 
to right): at the main sequence termination (TAMS), at the RGB tip, at the onset of thermal pulses on the AGB 
(Ist-TP), and at the last thermal pulse (AGB-tip). These phases comprehend the whole post-main sequence nuclear 
life of these stars. The bold solid lines (connected by dashed lines between stages), delimit the following boundaries 
(from top to bottom) between: the stellar surface and the mass lost by stellar winds, the H-rich envelope and the 
(H-exhausted) He core, and the He and CO cores. The chemical composition {X, Y) of the H-rich envelope is modified 
by the dredge-up (DU) episodes as marked. During the RGB and AGB phases, mass-loss significantly reduces the 
stellar mass (regions delimited by dots). See text for a complete discussion. 



ucts synthesized at He-shell flashes, mainly ^^C and ^^O 
of primary synthesis. 

In addition to the aforementioned dredge-up episodes, 
the surface chemical composition of TP-AGB stars with 
M[ <; 3.5M0 can be altered by the occurrence of 
hydrogen-burning via the CNO-cycle at the base of their 
convective envelopes (hot-bottom burning; HBB). This 
latter process also increases the surface abundance of ^He, 
and that of ^''N at the expenses of the newly dredged-up 
i^C (and possibly ^^O). 

Another crucial process characterising the evolution 
of low- and intermediate-mass stars is mass loss via stel- 
lar winds, both during the RGB (of low-mass stars) and 
AGB phases. From both theoretical and observational ar- 
guments we get indications that during the final stages of 
the AGB evolution of these stars, mass loss is able to al- 
most completely strip off their mantles. The final result is 
a white dwarf consisting almost completely of a C-0 core, 



and with very thin (in mass) superficial layers of H and 
He. 

It is just the combination of surface composition 
changes (due e.g. to dredge-up events) and mass loss that 
determines the contribution of these stars to the chem- 
ical enrichment of the interstellar medium (ISM). Such 
contribution is commonly quantified through the stellar 
yields, defined as the amount of newly synthesized ele- 
ments ejected into the ISM. 



3.2. Analytical derivation 

In order to guide our discussion. Figs. |l| and || schemati- 
cally present the basic structure of low- and intermediate- 
mass stars at some key-evolutionary stages. In our scheme, 
each star starts as a chemically homogeneous zero-age 
main sequence (ZAMS) configuration, and ends up as a 
CO white dwarf with virtually no H- or He-rich envelopes. 
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Intermediate-mass star 

.t. 




WD 



TAMS 



He-ignition 



Ist-TP 



AGB-tip 



evolutionary stage 

Fig. 2. The same as Fig. nlbut for an intermediate-mass star. In this case, the RGB phase is almost absent and with 
negUgible mass loss. 



Thus, the evolutionary stages depicted comprehend its 
whole post-MS nuclear evolution. 

Each stage is identified by the corresponding mass of 
the H-exhausted core Afc, namely: 

— the termination of the main sequence (TAMS). This 
is equivalent to the turn-off (TO) point of the corre- 
sponding isochrone (or SSP) with AfJ°; 

— the tip of the RGB (or the stage of central He ignition) 
with M<!^C!Bt. 

— the first thermal pulse on the AGB with M^"^^; and 

— the tip of the AGB with M^^^iBt (notice that M^'^^^ 
coincides with the remnant mass Afyvo)- 

Moreover, it is worth noticing that during the stages under 
consideration not only the mass-coordinate of the core(s), 
but also its chemical composition (either a He-core, or a 
CO-core surrounded by a very thin He-rich layer) can be 
easily singled out in stellar models. This can be appreci- 
ated looking at the lines depicted in Figs. [I] and g that 
schematically describe the evolution of the H-He and He- 
CO discontinuities. 

Outside the H-exhausted core, there is an envelope 
whose mass is reduced by both the core growth and mass 
loss, but can be temporarily increased (even if the former 
decrease is by far dominant) by dredge-up events. These 



latter events also change the envelope chemical composi- 
tion. We denote by Xj, Yj - with j = 0, 1, 2, 3 - the enve- 
lope fractional abundances (by mass) of hydrogen and he- 
lium respectively, where the superscripts indicate the orig- 
inal abundances if j = 0, and those after the j'^-dredge-up 
event if j > 0. 

In general, with the term "core mass" we denote the 
mass-coordinate of a specified chemical discontinuity (e.g. 
the mass of the H-exhausted core corresponds to the H-He 
discontinuity). This definition, however, cannot always be 
applied. For instance, stellar models at the TAMS do not 
show any sharp H-He discontinuity, but rather a chemical 
profile - left by either the radiative core or the recession 
of the convective core -, so that Mc does not corre- 
spond to any well-defined mass-coordinate inside the star. 
Anyhow, we can find a physically equivalent definition for 



MJ°, that is 



M. 



TO 



1 

-'^O JO 



Mi 



[Xo-X{Mr)]dMr, 



(6) 



where X{Mr) is the H abundance at the mass-coordinate 
Mr- In practice, Mc corresponds to the mean value of 
Mr across the chemical profile, weighted by the difference 
between the initial and local H abundance. We can also 
notice that the total mass of H burnt during the MS is 
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equal to Xq M^ ■ Moreover, Eq. (0) offers an alternative 
way to compute M^ from evolutionary tracks: 



M 



TO 



1 



^O^H 



LM,{t)dt 



(7) 



MS 



At this point we have defined (and illustrated in Figs. |] 
and|2|) all the ingredients necessary to estimate the nuclear 
fuel consumed in different evolutionary stages of low- and 
intermediate-mass stars. 

3.2.1. The total fuel 

Hereinafter we will briefly refer to total fuel as that con- 
sumed during the entire post-main sequence evolution. 
This can be evaluated with 



(Xo + 0.1)M, 



AGBt 



AGBt 



XqM. 



TO 



TO 



;iHc 



;iCO 



where (Xq + 0.1) M/'^'''- - Xq Mc^"^ represents the fuel 
required to convert a He-core of mass Mc into a CO- 
core of mass Mc^°^* (see Figs. and |); Ff"" and F^'^ 
represent the contributions to the fuel related to the stellar 
yields of ^He, and ^^C -I- ^^O elements (produced by He- 
burning reactions), respectively. In this paper, My{i) ¥\ 
denotes the stellar yield of element i (in mass units). 
The term referred to ^He is, by definition, simply 



Fy"'= = My(He) 
with 



My (He) 



[Y{t) 



-»'1^- 



(9) 



(10) 



where the integral is carried out over the entire stellar 
lifetime t*; Y{t) and dM/At are the current (at time t) 
surface ^He abundance and mass-loss rate, respectively. 

The newly synthesised CO nuclei from He-burning re- 
actions are clearly of primary origin. This applies, for in- 
stance, to ^^C and ^^O nuclei brought up to the surface 
by the third dredge-up, as they are synthesised at He- 
shell flashes during the TP-AGB evolution. Denoting with 
M^(CO) the corresponding primary yield (see Eq. p3| be- 
low), we can use the approximation: 



F^ 



CO 



= i.i(i-r')^'^y^(co) 



where 



Y' 



(1.1 



Yi. 



y')MP(co), 



X 



1,2 -t- -fl.2 



Yi. 



o.ir'AfP(co) (11) 



(12) 



This latter quantity involves H and ''He abundances at the 
onset of the TP-AGB phase, i.e. after the 2"'^ dredge-up 
or, if this latter does not occur as in lower mass stars, after 
the I''* dredge-up. Equation O accounts for the fact that 
a fraction Y' of Af^(CO) derives directly from nuclear 

^ For all the definitions related to yields, we refer to the 
classical work by Tinsley (1980), and to the recent calculations 
presented by Marigo (2001). 



burning of original helium, whereas the complementary 
fraction 1 — 1"' is synthesized starting from original hy- 
drogen through the sequence of both H- and He-burning. 

It is worth remarking that all the quantities that enter 
in the evaluation of Eq. (0) - i.e. the core masses Mj*-*, the 
remnant mass A/^^*^^* = AfwD, the chemical abundances 
after dredge-up episodes, Xj and Yj, and the stellar yields 
My (He) and M^(CO) - can either be easily derived from 
published stellar tracks, or are already tabulated in papers 
that deal with chemical yields. 

Finally, the two components of Eq. (H) can be also 
distinguished: 



AMh ~ Xo (M, 



AGBt 



A4^^)+My(He) 



+(i-y')*^v(co) 



(8) AM] 



Ho 



M, 



AGBt 



My^(CO) 



(13) 



(14) 



The above equations can be slightly modified when 
we consider stars (with A/j ^ S.SMq) experiencing hot- 
bottom burning (in addition to the dredge-up episodes) 
during the TP-AGB phase. In this case the dredged-up 
'^C (and possibly ^^O) may be converted into ^'^N, so that 
it is advisable to replace MP(CO) in Eqs. (|l|) - (|l|) with 
the total primary yield of the CNO elements, Af^(CNO). 

With the aid of Figs. H and 0, we can easily specify the 
contributions to the total fuel coming from the different 
post-main sequence shell(s)-burning phases: 



X^T — F-\ 



RGB 



F, 



CHoB+EAGB 



-Ftp- 



TP-AGB 



(15) 



where the rigt-hand side terms arc those derived in the 
following. 



3.2.2. The RGB fuel 

Frgb ^ (Ai+0.1)(A1 

-^(Fi - Fo) (Mi 



RGBt _ ^^ TON 



(16) 



A//, 



RGBt 



)• 



Here, the term related to the 1*** dredge-up is given by 
the amount of newly synthesized helium mixed up into 
the envelope mass M\ — M^ *, where M; is the initial 
stellar mass. We can notice that in this case the term 
Fy'~' = 0, as the first dredge-up involves material that has 
experienced only H-burning reactions. 

This equation requires the evaluation of Mc *, the 
core mass at He-ignition. This quantity is normally tab- 
ulated only for low-mass stars, because it critically de- 
termines the luminosity of HB and red clump stars. 
Moreover, in low-mass stars we always have Mc^*^^* — 
mJ''^ > 0, since the mass of the electron-degenerate core 
must grow up to 0.45 — 0.50 M© before He can ignite. 

For intermediate-mass stars, the RGB phase is practi- 
cally missing. Anyhow, for the sake of a uniform notation, 
we can still refer to i^RGB as the fuel (very little) consumed 
from the TAMS up to central He-ignition. Moreover, we 
notice that, in general, the quantity Mc * is not easily 
singled out in stellar tracks. In fact, in intermediate-mass 
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models the H-burning shell is rather thick (in mass), so 
that a sharp H-He discontinuity cannot be recognised. To 
overcome this difficulty we may adopt a simple approach, 
that is to derive M.^^^' from: 



M. 



RGBt 



M 



TO 



M, 



where 



Mdred ^ ^' ^° (Mi 



drcd 



M. 



TO 



)• 



(17) 



(18) 



In other words we assume that during the "RGB" phase 
the core is only affected by the l**' dredge-up, which re- 
duces its mass (by the amount Mdrod = Mc — Mc *), 
and increases the envelope helium content from Yq to Yi 
(normally a tabulated quantity). Actually the l'^' dredge- 
up usually proceeds inward across a chemical profile, in- 
stead of a H-He discontinuity implicitly assumed in the 
latter formula. Therefore, the above equations give just 
a crude approximation of Mc for intermediate-mass 

stars. 

3.2.3. The core He-burning -|- early-AGB fuel 

i^CHeB+EAGB ^ (^1 + 0.1) A'h^^^ - Xi M^"™' (19) 

+ {Y2 - Yl) (Afi - Afe'^P) . 



standard definition of stellar yields (and in Eq. 10). Then, 
the ^He contribution can be expressed as 



JV„ 



M, 



TP-AGB 



(He) = J2 [iYi+2 - 1^1,2) AMj 



(22) 



where Yj+2 corresponds to the abundance after the j 



th 



thermal pulse (or (j -I- 2)*"-dredge-up event), and AMj 
denotes the mass ejected during the j*^ pulse-cycle. The 
summation is performed over the total number, A'p, of 



pulse cycles. A functional form analogous to Eq. (22) 
should apply to M.J' (CO) as well, i.e.: 



Afo 



m;(CO)=^[xP(C 



0),+2 AM,] 



(23) 



where X^{C + 0)j+2 denotes the primary carbon and oxy- 
gen abundance in the envelope after the j^^ thermal pulse. 
It should be noticed that, in this case, the surface chem- 
ical composition at the onset of the TP-AGB phase does 
not contain any element of primary origin, i.e. the scaling 
term (corresponding to Yi 2 in Eq. E2h is set to zero in 
Eq. (g). 

Finally, we recall that in case of hot-bottom burning a 
better approximation is achieved using Ai?(CNO) instead 
of My (CO) (see earlier in this section). 



This expression can be easily understood when compared 
to Eq. ( 16|) . The term related to the 2"'^ dredge-up is zero 
whenever this does not occur (i.e. for M ^ 3 — 4 M©). 

Mc^ is normally a tabulated quantity, as in e.g. 
Marigo et al. (1998). 

Since core-He burning tends to convert the whole He- 
core into a CO-core, and since during most of the E-AGB 
phase the H-burning shell is off, one can also give a rough 
estimate to the fuel that comes from the E-AGB phase 
only: 



F, 



EAGB 



0.1 (Afc^^^ - A/e^''^*) 



(20) 



-(Y2-Yi)(Afi-A/, 



ITPn 



Finally, we remind the reader to consider the caution- 
ary remarks, expressed in Sect. 3.2.2, on the definition of 
Mc ' in the case of intermediate-mass stars. 



3.2.4. The TP-AGB fuel 



-Ftp- 



TP-AGB 



(Xi,2 + 0.1) [M 



AGBt 



M, 



1TP\ 



(21) 



-Af, 



TP-AGB 



(He) 



+(i.i-y')MP(co). 



It should be noticed that the quantity related to the 
stellar yield of helium, My'"^~^'^^(He), should be scaled 
with respect to the elemental abundance at the beginning 
of the TP-AGB phase, and not to the initial one as in the 



3.2.5. Further corrections 

The analytical relations presented so far consider the main 
contributions to the fuel, (i.e. nuclear burnings), but ne- 
glect other terms that possibly take part to the energy 
balance of a star, such as neutrino losses, and energy gains 
due to gravitational contraction. For the sake of complete- 
ness, these terms should be added to the right-hand side 
of Eq. (I2). However, neglecting them does not introduce 
a significant error in the evaluation of the fuel, as these 
terms only determine small corrections not exceeding a 
small percentage in most cases. 

Anyhow, we can easily get a rough estimate of the 
gravitational energy released by gravitational contraction 
during the evolution of low- and intermediate-mass stars 
with 



AEr 



_GMU 
RwD 



(24) 



which gives the order of magnitude of the gravitational 
energy of a white dwarf with a mass Mwd and radius 
^WD- This latter can be derived, for instance, adopting 
the mass-radius relation for white dwarfs (as derived from 
Tuchman et al. (1983) basing on Chandrasekhar (1939)): 



RwD = 0.019 (1 - 0.58 A/wd) 



(25) 



(expressed here in solar units) and setting AfwD = 



M. 



AGBt 



i.e. the mass of the core left at the end of the 



AGB. 
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Then, the corresponding contribution to the stellar fuel 
(in mass units) can be evaluated with 



AFa = - 



IAEg 
2 Ah 



(26) 



where the factor 1/2 accounts for the fact that only 
half of the gravitational energy goes in radiation (Virial 
Theorem) . A simple evaluation reveals that AFq is of or- 
der of just a few hundredths of Mq. The relative contri- 
bution of this term to the total fuel will be shown in a few 
examples discussed in Sect. |^. 

As far as neutrino losses are concerned, in the evolu- 
tionary phases under consideration they can be essentially 
of two kinds: 

1. Neutrino losses via (3 decays: they carry away a fraction 
of the energy produced by nuclear reactions. Since the 
constant An (Eq. g) includes only the energy that is 
effectively available to the star - thus excluding neutri- 
nos - no further correction is needed in our equations. 

2. Plasma neutrino losses: they cool the stellar cores dur- 
ing the RGB and AGB phases. However, the energy 
they carry away is always very small if compared to 
the stellar surface luminosity. In RGB stars, just some 
thousandths of the stellar energy leaves the star in the 
form of plasma neutrinos, whereas in the AGB this 
fraction is of a few hundredths- reaching a maximum 
of about 0.06 in the most luminous and massive AGB 
stars. Thus, they would lead to very small corrections 
- a few percent in extreme cases - to the fuel, that we 
will neglect in our approach. 

3.3. The case of massive stars 

Massive stars evolve through successive nuclear burn- 
ing stages in their interior up to the supernova (SN) 
explosion, possibly leaving either a neutron star or a 
black hole as remnant. In comparison to the case of low- 
and intermediate-mass stars, several factors hamper the 
derivation of simple and accurate formulas for the fuel 
consumed by massive stars, including the complex succes- 
sion of nuclear burnings in the latest stages of hydrostatic 
core evolution, and the complex energetics of the explosion 
event. 

Nonetheless, we may still provide simple formulas in 
order to give reasonable estimates to the fuel. First, we 
can suitably separate the total fuel into two components: 



^T = ^prc-SN + FsN 



(27) 



that refer to the pre-SN, and SN contributions, respec- 
tively. In analogy to Eq. (^, the first term can be gener- 
ahsed as: 



^(Xo + a;M)(Mc'-Afc'+') 



(28) 



i>2 



-XoAfc™+^«° + F,^°. 



The first term and the summation in the second term 
account for all the contributions to the fuel involved in 
the formation of the "onion-skin" structure inside the 
H-exhausted core. Here the increase of the index i cor- 
responds to more and more advanced nuclear burnings 
(e.g. 1 = 1 refers to H-burning, ? = 2 is related to He- 
burning, etc.) and M^ denotes the mass-coordinate of the 



,-th 



-element exhausted core (e.g. M^ ° for z = 1, M, 



CO 



for i = 2, etc). In general, the quantity Xij is related to 
the amount of energy generated in the conversion of nuclei 
i into nuclei j, normalised to that provided by the conver- 
sion of H into "^He. It can be estimated from the difference 
between the binding energies per nucleon, Ba, that char- 
acterize the nuclei involved (i = final, j = initial): 



B 



A(t) 



B 



AU) 



B. 



He 



(29) 



with respect to that of ^He (the main product of H- 
burning). We can notice that, by definition, Bn = 0, 
and that a;2,i — {Bi2q — Bai^^) / Bi^^ '^ 0.1, as already 
seen in the previous sections. In other words, Xi^i repre- 
sents the global efficiency - relative to that of H-burning 
- of all nuclear burnings (from that of He to the i*'^-one) 
which have successively taken place in the shell of mass 



(Me 



M. 



i+1 



TO 



Finally, the term -XqMc^"^ + F^'' + F^° has the 
same meaning as in Eq. (||), but with a difference: here 
the chemical yields refer only to those produced during 
the pre-SN evolution, that is, to the so-called wind yields. 

The practical applicability of Eq. (|2q) is determined by 



the fact that the quantities Mc °, Mc , and the wind 
yields, are normally presented in published pre-SN models 
(see e.g. Maeder 1992; Woosley & Weaver 1995; Portinari 
et al. 1998). Conversely, the interior "onion-skin" struc- 
ture is generally not explicitly given. Anyway, we can use 
simple approximations to derive lower and upper limits 
to the second right-hand side term of Eq. (|2^). A lower 
limit is obtained neglecting any nuclear conversion into 
elements heavier than CO nuclei (a;2,i = 0.1), so that the 
term simply reduces to (A"o-l-0.1)Mc . An extreme upper 
limit corresponds to supposing that the He-exhausted core 
has been wholly converted into the iron-group elements. 
In this case, we have Xi^i = x(^^C — >^^ Fe) ~ 0.16 and the 
second right-hand side term becomes {Xq -\- 0.26)Mc 

Therefore, these simple estimates suggest that -Fprc-SN 
should be comprised within a rather narrow range, of 
width ~ 0.16Afc . This narrowness is determined by the 
relatively low amount of energy that is available from the 
conversion of CO into heavier elements. 

Computing _Fsn can be far more difficult, because of 
the complicated energetics involved in the SN explosion. 
Anyway, a simple order-of-magnitude estimate comes from 
the conversion of the typical energy irradiated by SNe, 
XO^^ erg, into the equivalent mass of burnt H: 

FsN~ 0.1 Mq. (30) 

Therefore, we would expect that usually i^sN *C -fprc-SN- 
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Fig. 3. Accuracy checks. Stellar fuel consumption as a function of the stellar mass, calculated for two sets of stellar 
models with different initial metallicity (as indicated) and considering various post-MS phases (labelled nearby the 
corresponding curve) . The results predicted by the "exact" method given by Eq. (||) is compared with those obtained by 
adopting the analytical formalism here proposed (dashed lines) . The total fuel is calculated both without (long-dashed 
line) and with (short-dashed line) the gravitational term expressed by Eq. (26). 



4. Illustrative applications 

Here we will apply the analytical formalism derived in the 
previous section to a few representative examples, in the 
domain of low- and intermediate-mass stars. First, we will 
check the accuracy of the formulas, applying them to a 
particular case in which we know in advance that the stel- 



lar inputs are homogeneous and self-consistent (Sect. 4.1). 
For this specific case, we will analyse the fraction of the 



fuel that is in the form of chemical yields (Sect. 4.2). This 
will also allow us to discuss the general features of the 
fuel as a function of the stellar mass and metallicity, and 
to quantify the error we may introduce if we neglect the 
effect of dredge-up episodes (and, in general, of chemical 
yields) in estimates of the fuel consumption. Finally, we 
will illustrate a case in which heterogeneous stellar data 
are used to model the emitted light and chemical yields 
(Sect. p~^ ). This will allow us to remark on the consis- 
tency requirements for the input stellar data in chemical- 
spectrophotometric models of stellar populations. 



4.1. Accuracy checks 

The accuracy of the analytical formalism presented in 
Sect. H should be tested by comparing its predictions 
to the corresponding results obtained with the "exact" 
method already mentioned in Sect, g, which is essentially 
based on the calculation of the integral in Eq. (|^). Of 
course, this can be done only if we have a homogeneous 
set of stellar models, providing tables with both (i) the 



evolutionary tracks (or isochrones), and (ii) the chemical 
yields and core masses. 

Here we adopt a set of low- and intermediate-mass evo- 
lutionary models with 0.9Mq ^ M < 5Mq, Zq = 0.008, 
Yq — 0.25 computed by Girardi et al. (2000) and Marigo 
et al. (1999), which follow the evolution from the ZAMS 
up to the end of the AGB, and predict the correspond- 
ing stellar yields (Marigo 2001). These models are more 
extensively described in Appendix A. 

The outcome of the test is shown in Fig. ||, for two 
different initial metallicities {Zq — 0.008 and Zq = 0.019). 
It turns out that the analytical prescriptions reproduce 
remarkably well the fuel calculated with the "exact" 
method, the relative difference ranging typically within 
5 %, mostly concerning the RGB, TP-AGB, and the to- 
tal post-MS contributions. As expected, the largest dif- 
ferences result in the analytical derivation of -Frgb and 
^CHcB+EAGB of intermcdiatc-mass stars, due to the dif- 
ficulty of defining the core mass-coordinate Mc * (see 
remarks in Sect. 3.2.2). 

Figure S displays also the effect of the correction term, 
Fq, calculated with Eq. (Eq), which is added to the total 
fuel derived from Eq. ^. As expected, the convergency 
towards the "exact" result is improved by a small amount 
(few percentiles at maximum) with the inclusion of the 
gravitational term. Finally, we notice that even if some 
of the partial terms (i.e. Frgb, ^HcB+eagb) may be not 
well reproduced, the analytical prescription for the total 
post-MS fuel Ft is quite accurate. This reflects the small 
number of stellar parameters involved in Eq. (ph . 
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Fig. 4. Upper panel: Total post-main sequence fuel con- 
sumption by low- and intermediate-mass models with ini- 
tial composition (Yq = 0.25, Zq — 0.008), calculated with 
and without the terms iF^'^ and FS''^) related to the 
chemical yields in Eq. (H). Middle panel: masses of the 
H-exhausted cores as defined in the text. Bottom panel: 
chemical yields of ''He and primary CNO. See text for 
more details 



4.2. The effect of surface chemical changes 

Figure |j displays the predicted stellar fuel consumed over 
the post-main sequence evolution, as a function of the 
initial stellar mass, again for the Girardi et al. (2000) 
and Marigo et al. (1999) tracks. For illustrative purposes, 



we first compute the fuel according to Eq. (m but ne- 
glecting the terms due to the dredge-ups (i.e. setting 
FyH'= = FyC'O = 0). This is referred to as case A. The in- 
clusion of these "chemical" terms corresponds instead to 
case B, that is the correct derivation of the fuel. 

As illustrated in Fig. H the differences between the two 
cases are remarkable for this particular set of stellar mod- 
els. As expected, the total fuel consumption according to 
case A is systematically underestimated with respect to 
case B. In particular, the following points are worthy of 
notice. 

— The discrepancy - between the results obtained in case 
A and the correct ones with prescription B - increases, 
on average, with the stellar mass, reflecting the trend 
of the chemical yields (bottom panel). 

— Lower mass stars (with A/ ~ 1 Mq) experience only 
the first dredge-up during their evolution, so that the 
"chemical" terms almost do not contribute to their Ft . 

— At around M ~ SM© both the actual total fuel Ft - 
calculated according to case B -, and the yield of he- 
lium, A/y(He), present a maximum, which corresponds 
to the longest duration of the TP-AGB phase for this 
set of stellar models (see Marigo 2001). 

— At larger stellar masses, prescription A predicts a dras- 
tic drop of the fuel, which does not show up with 
prescription B. This can be explained comparing the 
curves of MJ^ and Afc"^^^* (middle panel of Fig. |). 
The difference (Afc * — Afc ) progressively reduces 
at increasing stellar mass, eventually becoming nega- 
tive. In other words, in the most massive models the 
final cores left at the end of the AGB phase may have 
even lower masses that the cores built up at the end 
of the MS. In these cases, the fuel calculated with pre- 
scription A would eventually attain negative values, 
thus losing significance. This example is already a clear 
indication that the evaluation of the stellar fuel sim- 
ply basing on core masses is not a correct procedure, as 
it may even produce meaningless results. Conversely, 
looking at predictions obtained with prescription B, 
we can see the total fuel first decreases and then flat- 
tens towards the most massive models (keeping always 
positive), due to the combined effect of both Mc and 
chemical yields. 

The weight of the "chemical terms" can be appreci- 
ated also in Fig. ^ which displays the results of the ana- 
lytical estimation of the fuel consumed during the RGB, 
TP-AGB and total post-MS phases, as a function of the 
initial stellar mass and metallicity. We can see that the 
effect of the I*'' dredge-up occurring on the RGB is rather 
small, whereas accounting or not for the surface chemical 
changes during the AGB leads to quite different results, 
i.e. the fuel can be substantially underestimated if not 
doing correctly. Moreover, we notice that the stellar fuel 
consumption presents a marked metallicity dependence, 
i.e. it increases, on average, at decreasing metallicity for 
all the evolutionary phases here considered. 
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Fig. 5. Stellar fuel consumption (derived with the analytical formalism) of low- and intermediate-mass models with 
different choices of the initial metallicity, and referring to i) the RGB phase (left panels); TP-AGB phase (middle 
panels); and total post-MS evolution (right panels). Results of calculations including the effect of dredge- up events 
and hot-bottom burning (solid lines) are compared to the case in which such effect is neglected (dotted lines). See text 
for more details 



It is also clear that the stellar fuel - calculated with 
the proposed formalism - is positively correlated to the ef- 
ficiency of the dredge-up episodes. To give an example, let 
us consider the TP-AGB evolution of a given stellar model. 
For the sake of simplicity, we assume that at each ther- 
mal pulse a dredge-up occurs with a constant efhciency, 
A = AMdrcd/AA/c, defined as the fraction of the core mass 
increment during an inter-pulse period which is dredged- 
up into the envelope at the subsequent thermal pulse. This 



implies that, every time a dredge-up takes place, the core 
mass is reduced by the amount AAA/c. The TP-AGB fuel 
consumption can be consequently expressed as: 

(],,f AGBt _ 11^ ITP^, 

Ftp„agb ^ (^1,2 + 0.1) ^'''^ , ''''' 1 (31) 

u - -^j 

that is equivalent to Eq. (pl|). 

It follows that, if we neglect the effect of third dredge- 
up, the stellar TP-AGB fuel consumption would be un- 
derestimated by a factor 1/(1 — A). Setting A = 0.50, as 



12 



Marigo & Girardi: Coupling emitted light and chemical yields from stars 



adopted in the calculations presented in Fig. H, the un- 
derestimation factor is ~ 2, i.e. the actual TP-AGB fuel 
would be twice larger than that one would obtain setting 
A = 0. Finally, it is worth noticing that such effect may be 
even larger if one assumes the extreme dredge-up efficien- 
cies (A ^ 1) of recent TP-AGB calculations (i.e. Herwig 
2000). 

Of course, only case B represents the correct ap- 
proach to be followed in order to estimate the stellar fuel. 
However, discussing case A is still useful. In fact, case A 
represents a sort of lower limit to the fuel, that depends 
only on the core masses at the MS and AGB termination 
stages. Since Mc is a monotonic function of the initial 
mass Afj, and Mc ' — Mwd, it follows that this lower 
limit is intimately linked to the MwD(Afi) relation (see 
also Girardi & Bertelli 1998). This latter is known as the 
initial-final mass relation (IFMR) , and is well constrained 
by observations of nearby white dwarfs (both in the fields 
and in open clusters). Thus, the observed IFMR poses a 
direct constraint ~ i.e. a lower limit - to the behaviour of 
the stellar fuel as a function of the mass. 

In the context of galaxy models, case A can be seen 
as the fuel (intended as the luminosity integral of Eq. ||) 
that we should be dealing with in our spectro-photometric 
model - instead of the one given by case B - under the 
assumption that low- and intermediate-mass stars do not 
contribute to the chemical evolution of the system under 
consideration. Needless to say, this situation is not realistic 
at all. 



4.3. The need for a consistency choice in galaxy 
models 

The analysis and relative discussion carried out in the pre- 
vious sections should have already convinced the reader on 
the intimate connection between light emission and chem- 
ical yields from stars. Up to now, we have illustrated this 
fact using a homogeneous set of stellar models (Girardi et 
al. 2000; Marigo et al. 1999, and Marigo 2001). To make 
concepts even clearer, let us suppose to opt for an incon- 
sistent choice in modelling a galaxy, that is combining the 
TP-AGB tracks of Marigo (2001, M2K; see Appendix A) 
to model the spectrophotometric evolution, with the stel- 
lar yields of Renzini & Voh (1981, RV81) to model the 
chemical evolutionn 

In this case, we can get two different estimates of the 
fuel: the first one is derived from our analytical formulas 
applied to RV81 yields and core masses and corresponds 
to fuel adopted in the chemical evolution model, and the 
second one comes from the exact method applied to M2K 
evolutionary tracks and corresponds to fuel adopted in the 
spectrophotometric model. Our primary aim is to see how 



* Together with those of van den Hoek & Groenewegen 
(1997), M2K and RV81 yields are commonly used in chemi- 
cal evolution models of galaxies. A detailed comparison among 
the different set of stellar yields is provided in Marigo (2001), 
to whom we refer for all details. 



they compare. If the two different inputs were consistent, 
the two estimates of the fuel should agree to within a 
few percent, just as in the case of Fig. 0. Otherwise, this 
exercise would allow us to quantify the degree of mismatch 
(inconsistency) between these two sources of stellar data. 

We use Equation (|2|) to derive the fuel consumption 
of TP-AGB models with initial metallicity Zq = 0.02 pre- 
sented by RV81. The results are shown in Fig. |q (left 
panel). We adopt the final core masses (Mc ) and 

chemical yields presented in tables 3f and 3e of RV81, 
this latter table corresponding to calculations of TP-AGB 
models (with M > 3.3 AIq) with hot-bottom burning car- 
ried out with the mixing-length parameter a — 2.0. For 
RV81 models, the amounts of mass lost during the RGB 
phase of low-mass stars (to be properly subtracted from 
the total yields to get the TP-AGB yields), the envelope 
abundances of ^He (necessary to calculate the quantities 
Y'), and the corresponding core masses at the onset of the 
TP-AGB phase (Mc^ ) are derived following the recipes 
described in their section 2. 

For M2K models with Zq = 0.019 and a = 1.68, the 
TP-AGB fuel has been already presented in our previous 
Sects. 11 (see Fig. ^. The left panel of Fig. |^ shows, 
instead, M2K fuels corresponding to the same value of 
a = 2.0 as adopted in the RV81 models presented in the 
same figure. 

For both sets of stellar models, we also compute the 
TP-AGB integrated bolometric luminosity as a function 
of SSP age, according to: 



L 



TP-AGB 
bol 



^ An b{t) Ftp^agb{Mto) 



(32) 



(where all quantities are defined in Sect. 0). They are dis- 
played in the right panel of Fig. ||. 

It is clear, from the comparison between RV81 and 
M2K results, that the TP-AGB fuel of RV81 models is 
much larger than the one of M2K models, for masses Mi ^ 
3Mq, corresponding to stellar ages shorter than about 
300 Myr. More modest but still non-negligible differences 
are also present at lower masses (older ages). 

In addition to the TP-AGB fuel as a function of the 
initial stellar mass. Fig. y (left panel) shows also the con- 
tribution of the "chemical" terms to Ftp-agb (i-e. the 
second and third terms in Eq. ^l]) for both sets of models. 
It turns out that i) in both cases the yield contribution to 
the fuel is generally relevant (except for the lowest mass 
models that do not experience both the third dredge-up 
and hot-bottom burning), ii) the differences between RV81 
and M2K are still considerable. We can notice that for 
Mi > 4 M0 the yield contribution to the TP-AGB fuel in 
RV81 is even larger than the total TP-AGB fuel in M2K 
models. 

As expected, the trend of -Ftp-agb as a function 
of Mi (left panel of Fig. y) reflects in the behaviour of 
^boi~ ^^ ^ function of the corresponding SSP age (for 
Mtq — M\] right panel of Fig. ||). In fact, here the most 
relevant differences between RV81 and M2K predictions 
show up at the youngest ages (hence higher Mxo)j i-e. 
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Fig. 6. Left panel: TP-AGB fuel calculated with the analytical formalism here presented and applied to stellar models 
of Renzini & Voli (1981, RV81; dashed line), and Marigo (2001, M2K; continuous line). Right panel: Contribution 
(in arbitrary units) of the TP-AGB phase to the integrated bolomctric luminosity of SSPs as a function of age. For 
RV81 models the evolutionary flux b{t) is calculated following Maraston (1998), whereas for M2K models we base on 
Girardi et al. (2000) evolutionary calculations, in both cases with the assumption of the Salpeter IMF. See text for 
more details. 



the luminosity contribution of the TP-AGB phase accord- 
ing to RV81 models is much greater than predicted by 
M2K models. 

On the basis of the above comparative analysis, let 
us now consider what happens if M2K tracks are used to 
model the spectro-photometric evolution, whereas RV81 
yields are used to model the chemical evolution. The most 
striking inconsistency in this model would show up at 
Mi ^ AMq (or equivalently at logi ^ 8.4), where the 
fuel necessary to account for RV81 yields would already 
be larger than the fuel necessary to account for the lumi- 
nosity of M2K models. For TP-AGB stars in such a range 
of masses and ages, our galaxy model we would be either 
strongly overestimating the chemical yields or, equiva- 
lently, significantly underestimating the emitted light. The 
degree of inconsistency is smaller, but still significant, at 
smaller masses (longer ages). 

Since the TP-AGB phase provides a significant fraction 
of the total energy emitted by stellar populations (meant 
either as bolometric hght, or as chemical yields), one may 
deduce that such models can hardly be considered as self- 
consistent in terms of their energetics. 

Moreover, a number of other subtle inconsistencies 
may affect such a galaxy model. For instance, RV81 and 
M2K models predict quite different initial-final mass re- 
lations (see M2K for a discussion). Whereas a chemical 
model that uses RV81 yields is expected to produce a 
relatively large fraction of "massive" white dwarfs (with 
-^WD ~ l-^©), a spectro-photometric model that adopts 



M2K tracks, under the same assumptions for star forma- 
tion rate and age-metallicity relation, is expected to pre- 
dict much fewer of them. Then, in view of modelling, for 
instance, the formation of galactic halos one would face 
the embarrassing question: which one of the two theoret- 
ical white dwarf mass distributions should be considered 
for a comparison with the observational data of field WDs 
and MACHOs? Any is the answer, it will be the wrong 
one. 



5. Final remarks 

This study has highlighted the intimate relation between 
the stellar emitted light and chemical yields. It has been 
clarified that a significant fraction of fuel consumption 
(emitted light) may be eventually "deposited" in the form 
of He and CO yields. This fraction can be quantified by 
means of simple analytical relations. We have illustrated 
this in detail for low- and intermediate-mass stars, also 
suggesting how the same analysis could be extended to 
massive stars. Then, the derived formalism offers a useful 
consistency check that can be applied to galaxy models 
whenever different sources of stellar data are employed to 
model the chemical and spectrophotometric evolution. 

Such consistency check expresses essentially the basic 
condition of energy conservation, which can be read as: 
the emitted light implies a precise amount of stellar nucle- 
osynthesis, part of which corresponds to chemical yields. 
Other fundamental aspects, such as mass conservation. 
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are already included in the basic equations of chemical 
evolution models. 

It is important to remark that the present analysis 
is general, and independent of the specific set of stellar 
tracks/yields in use. As a matter of fact, analysing the 
quality of any given set of stellar (or galaxy) models, is 
beyond the scope of this paper. 

Anyway, a few remarks can be made. As indicated by 
our analysis, the primary yields of He and CNO elements 
are those more directly related to the emitted light of a 
stellar population. In galactic chemical evolution models. 
He and CNO yields critically determine the evolution of 
Y{t) (or, equivalently, the AY/AZ ratio) and Z{t) rela- 
tions. Our formalism suggests that, if a given set of stellar 
models predicts a particular evolution of the integrated 
light, the same set of models necessarily implies a partic- 
ular evolution of Y{t) and Z{t). 

Of course, the situation is somewhat complicated by 
the fact that most of the metals are produced by massive 
stars, which are characterised by rapid evolution, so that 
soon they do not contribute to the integrated light any- 
more. Hence, light and chemical evolution may be dom- 
inated by different stars, possibly justifying the approx- 
imation that photometric and chemical evolution could 
be treated separately. Actually, this assumption may be 
safely applied to some cases (e.g. in the study of the 
present Galaxy), but becomes questionable (and risky) 
once galaxy m,odels are constructed to describe the evolu- 
tion of galaxies over most of their history, since their for- 
mation up to present times. An example of this "unifying" 
approach is given by the studies of high-redshift galax- 
ies based on models "calibrated" on local galaxy samples. 
Clearly, our analysis indicates that the energetics neces- 
sary to explain the light of the distant sample has strict 
implications for the chemistry now observed in the local 
sample, and vice-versa. 

Finally, we would like to conclude outlining the main 
points of the present work: 

— The stellar fuel (emitted light) consumed by a star 
during a specified evolutionary phase may be schemat- 
ically decomposed as the sum of two terms: a contri- 
bution locked in the star and related to the mass of 
the core regions and envelope composition, and a con- 
tribution emitted by the star during that phase in the 
form of chemical yields. These terms are explicited in 
this paper by means of suitable and simple analytical 
formulas, covering various phases of stellar evolution. 

— When coupling the spectro-photometric and chemical 
evolution of a system, we should care that the basic 
contributions of the stellar component - i.e. light and 
newly synthesised elements - are included in the mod- 
els consistently one with each other, preferably derived 
from the same set of stellar models. If this is not the 
case, i.e. the input stellar data come from heteroge- 
neous sources, we could at least measure the degree 
of internal mismatch of the model. The analytical for- 
malism developed in this work can serve the purpose. 
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Appendix A: Isochrones with improved TP-AGB 
models 

Girardi et al. (2000) presented a large set of evolution- 
ary tracks for low- and intermediate-mass stars, covering 
the mass range from 0.15 to 7 Mq, for six different initial 
chemical compositions. These models account for a moder- 
ate amount of convective overshooting from stellar cores, 
which makes the upper mass limit for the development 
of the AGB phase to be about 5 Mq , instead of 8 Mq as 
predicted by classical stellar models. 

Since the evolutionary calculations, carried out by 
means of the complete stellar code, stopped after the first 
few thermal pulses, Girardi et al. (2000) followed the sub- 
sequent TP-AGB phase with the aid of a simple synthetic 
algorithm (cf. Girardi & Bertelli 1998). In this way, sets 
of isochrones including the complete TP-AGB phase were 
presented. 

That approach to the TP-AGB evolution, however, is 
far too simple if compared to the detailed semi-analytical 
TP-AGB models presented by Marigo et al. (1999, and 
references therein). These latter models include the most 
crucial aspects of TP-AGB evolution, such as the third 
dredge-up, hot-bottom burning, deviations from the core 
mass-luminosity relation, and mass-loss rates related to 
the pulsational period. Moreover, the basic model param- 
eters have been calibrated such as to reproduce the carbon 
star luminosity functions in the LMC and SMC. Chemical 
yields from these models are presented by Marigo (2001). 

Marigo et al. (1999) tracks start at the first ther- 
mal pulse computed by Girardi et al. (2000), so that 
both sets of models are perfectly complementary. They 
follow the entire TP-AGB evolution up to the complete 
ejection of the stellar envelope. These synthetic TP-AGB 
tracks are available, so far, only for chemical compositions 
[Zo = 0.004, Fo = 0.240], [Zq = 0.008, Yq = 0.250] and 
[Zo = 0.019, yo = 0.273]. For these metallicities, the syn- 
thetic TP-AGB models have been suitably combined to 
the tracks calculated by Girardi et al. (2000), so that the- 
oretical isochrones have been constructed using the same 
method as in Girardi et al. (2000). 

Thanks to their continuity and homogeneity, the 
present isochrones are well suited for the population syn- 
thesis of galaxies. Moreover, they perfectly meet the con- 
sistency requirements discussed in this paper: the emitted 
light computed with the present isochrones (with the aid 
of Eq. Q) , would be completely consistent with the chem- 
ical evolution derived from the tables of chemical yields 
and remnant masses in Marigo (2001). 

Moreover, it is worth emphasizing that Marigo et al. 
(1999) models reproduce a series of observational proper- 
ties of AGB stars in nearby galaxies, and of nearby white 
dwarfs (see Marigo 2001 for details). Therefore, one may 
expect that they are adequate to model the integrated 
light emitted by TP-AGB stars. 
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Fig. A.l. The fractional contribution of luminous AGB 
stars tho the integrated light of LMC and SMC star 
clusters, as inferred from Frogel et al. (1990) data (full 
dots with error bars), as a function of the equivalent 
SWB types. This latter represents a rough age ranking 
for LMC and SMC clusters, that corresponds to the log- 
arithmic age scale presented in the upper axe. The same 
quantity is plotted, as a function of age, as derived from 
our isochrones of metallicity Zq — 0.004 (dashed line), 
Zo = 0.008 (continuous fine) and Zq = 0.019 (dotted line). 



A.l. Integrated ligth contribution from AGB stars 

In order to check this latter point, we compare the present 
models to the data from Frogel et al. (1990), who measured 
the contribution of luminous AGB stars to the integrated 
light of Magellanic Cloud star clusters of various ages. To 
make the comparison, we proceed as follow: First, from 
their table 1 we compute the bolometric light coming from 
stars with il/boi < —3.8 in each cluster (assuming distance 
moduh of 18.5 and 18.9 mag for the LMC and SMC, re- 
spectively). This Mboi limit guarantees that we are dealing 
only with AGB stars brighter than the RGB-tip. Second, 
we divide the clusters into bins of equivalent SWB (Searle 
et al. 1980) types (cf. table 4 in Frogel et al.). Third, 
we compute the ratio, £boi(.^^boi < — 3.8)/Lboi(total), be- 
tween the luminosity coming from stars with Afboi < —3.8 
and that emitted from the whole cluster (cf. their table 4) . 
Finally, we estimate the corresponding error bars from the 
total number of Mboi < —3.8 stars in each SWB bin. 



The results are presented in Fig. |A.1| , which also dis- 
plays the approximate relation between SWB type and the 
logarithm of cluster age, i.e. log(i/yr) ~ 6.66 -I- 0.48 SWB. 
The observed data in this figure are indeed very sim- 
ilar to the original ones in figure 15 of Frogel et al. 
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(1990). Superimposed on the data, we plot our theo- 
retical expectation for the same quantity, iboi(Afboi < 
— 3.8)/iboi(total), as derived from the present isochrones. 

From the figure, it is well evident that the models cor- 
rectly reproduce the general behaviour of the relative lu- 
minosity contribution over the complete age interval in 
which AGB stars are present (i.e. 8 ^ log(i/yr) <^ 10). In 
particular, the minimum values at log(i/yr) ~ 8 and ~ 10 
are clearly present in the models, as well as the maxi- 
mum at log(i/yr) ~ 8.6, that would correspond to the 
longest TP-AGB hfetimes (see Marigo 2001). It is also 
worth noticing that according to the models at decreasing 
metallicity the luminosity contribution from AGB stars 
increases. 

For ages younger than 10^ yr, we have computed the 
same quantities using Bertelli et al. (1994) isochrones. The 
integration of LboK-^boi < —3.8) was limited to stars with 
logTcff < 3.8, in order to exclude the blue main sequence 
stars from this term. As shown in Fig. A.l, the models 
predict values Lboi(A^boi < — 3.8)/Lboi(total) of around 
0.1 at these young ages. They are in reasonable agreement 
with observations, despite of the fact that the empirical 
data present a slight hint for an increase in -Lboi(-^boi < 
—3.8) at even younger ages (^ 10^ yr) (we notice, however, 
this feature is based on too few stars to be completely 
rehable). 

With the cautionary remarks that (i) the SWB clas- 
sification represents just a crude age ranking for clusters, 
whose consistency with our models has not been checked 
for, and (ii) that there is plenty of room for improving the 
empirical data for AGB stars in clusters, we conclude that 
our models reproduce the observed contribution of AGB 
stars to the integrated light of SSPs in a satisfactory way. 

A.l. Useful data 

Complete tables with the improved isochrones de- 
scribed in the previous section are available upon 
request to the authors, and at the WWW site 
http: //pleiadi .pd. astro . it. 

For the present models, Tables Al - A3 contain the 
stellar fuels as a function of the initial mass, for various 
post-MS phases. All fuels have been calculated by means 
ofEq. M. 

We also present fitting relations to Mc as a function 
of the initial mass M; and metallicity Zq, derived from the 
stellar models here considered (with overshooting) : 



Table Al. Post-MS stellar fuels (Mq) - Zq ^ 0.019 



M 



M. 



TO 



TO 



0.463 - 0.654 M; -I- 0.338 M^ (A.l) 

-0.024 log(Zo/0.019) for M; <, 1.4 
0.006 -t- 0.115 Afi -I- 0.023 M^ (A.2) 

-0.057 log(Zo/0.019) for Afi > 1.4, 



and from classical models (without overshooting) with ini- 
tial solar composition {Yq = 0.273 , Zq = 0.019): 



Mi 


^RGB 


-FcHob+E-AGB 


f^TP-AGB 


Ft 


0.868 


2.371E-01 


l.OOlE-01 


2.892E-02 


3.661E-01 


0.934 


2.387E-01 


1.014E-01 


1.400E-02 


3.542E-01 


1.005 


2.405E-01 


1.023E-01 


1.912E-02 


3.620E-01 


1.082 


2.428E-01 


1.059E-01 


1.887E-02 


3.675E-01 


1.163 


2.385E-01 


1.042E-01 


2.402E-02 


3.667E-01 


1.248 


2.287E-01 


1.061E-01 


2.778E-02 


3.625E-01 


1.334 


2.141E-01 


1.091E-01 


3.049E-02 


3.537E-01 


1.420 


1.951E-01 


1.095E-01 


3.577E-02 


3.405E-01 


1.504 


1.731E-01 


1.158E-01 


4.829E-02 


3.372E-01 


1.588 


1.573E-01 


1.182E-01 


5.916E-02 


3.347E-01 


1.672 


1.415E-01 


1.179E-01 


7.404E-02 


3.334E-01 


1.756 


1.230E-01 


1.245E-01 


8.503E-02 


3.326E-01 


1.839 


1.006E-01 


1.327E-01 


1.065E-01 


3.398E-01 


1.923 


6.746E-02 


1.483E-01 


1.222E-01 


3.380E-01 


2.000 


1.910E-02 


1.870E-01 


1.485E-01 


3.546E-01 


2.200 


1.640E-02 


1.781E-01 


1.766E-01 


3.711E-01 


2.500 


1.130E-02 


1.652E-01 


2.145E-01 


3.909E-01 


3.000 


8.800E-03 


1.597E-01 


2.693E-01 


4.378E-01 


3.500 


8.600E-03 


1.604E-01 


2.149E-01 


3.839E-01 


4.000 


9.700E-03 


1.780E-01 


1.363E-01 


3.239E-01 


4.500 


1.090E-02 


1.944E-01 


l.OOlE-01 


3.052E-01 


5.000 


1.170E-02 


2.144E-01 


1.030E-01 


3.289E-01 


Table A2. Post-MS stellar fuels {Mq) - Zq = 0.008 


Mi 


i^RGB 


-fcHeb+E-AGB 


i^TP-AGB 


Ft 


0.850 


2.542E-01 


9.550E-02 


4.851E-02 


3.982E-01 


0.918 


2.522E-01 


9.746E-02 


4.888E-02 


3.986E-01 


0.992 


2.516E-01 


1.006E-01 


3.666E-02 


3.889E-01 


1.071 


2.514E-01 


1.017E-01 


4.052E-02 


3.936E-01 


1.154 


2.499E-01 


1.048E-01 


5.144E-02 


4.061E-01 


1.239 


2.414E-01 


1.066E-01 


5.500E-02 


4.030E-01 


1.407 


2.061E-01 


1.055E-01 


7.898E-02 


3.906E-01 


1.499 


1.800E-01 


l.lOlE-01 


9.679E-02 


3.867E-01 


1.583 


1.626E-01 


1.138E-01 


1.172E-01 


3.936E-01 


1.667 


1.452E-01 


1.168E-01 


1.392E-01 


4.012E-01 


1.750 


1.235E-01 


1.224E-01 


1.644E-01 


4.102E-01 


1.832 


9.532E-02 


1.347E-01 


1.985E-01 


4.285E-01 


1.900 


2.450E-02 


1.866E-01 


2.317E-01 


4.427E-01 


2.000 


2.380E-02 


1.850E-01 


2.451E-01 


4.540E-01 


2.200 


1.720E-02 


1.808E-01 


2.758E-01 


4.738E-01 


2.500 


1.260E-02 


1.759E-01 


3.047E-01 


4.933E-01 


3.000 


1.020E-02 


1.871E-01 


3.539E-01 


5.513E-01 


3.500 


9.700E-03 


2.116E-01 


2.581E-01 


4.794E-01 


4.000 


1.030E-02 


2.288E-01 


1.552E-01 


3.943E-01 


4.500 


1.080E-02 


2.426E-01 


1.354E-01 


3.888E-01 


5.000 


1.150E-02 


2.589E-01 


1.083E-01 


3.787E-01 








for Afi 


<1.4 


Afe^° 


~ 0.015 + 0.063 Afi + 0.019 Aff 


(A.4) 








for Afi 


> 1.4. 



All masses are expressed in solar units. We remind that 
according to our definition of Afc (see Eq. M), the fuel 
burnt during the MS can be estimated from: 



M 



TO 



0.123 + 0.066 Mi - 0.037 Af;^ 



(A.3) 



Fms — Xq M, 



TO 



(A.5) 
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Table A3. Post-MS stellar fuels (M©) - Zo = 0.004 



Mi -Frgb -FcHob+E-AGB -Ftp-agb Ft 

0.817 2.536E-01 9.664E-02 5.876E-02 4.089E-01 

0.888 2.527E-01 1.005E-01 6.502E-02 4.182E-01 

0.966 2.517E-01 1.040E-01 6.470E-02 4.204E-01 

1.046 2.499E-01 1.067E-01 7.085E-02 4.274E-01 

1.131 2.467E-01 1.091E-01 8.151E-02 4.373E-01 

1.218 2.397E-01 1.115E-01 9.972E-02 4.510E-01 

1.307 2.223E-01 1.122E-01 1.392E-01 4.736E-01 

1.396 2.002E-01 1.147E-01 1.771E-01 4.919E-01 

1.481 1.741E-01 1.153E-01 1.989E-01 4.882E-01 

1.564 1.508E-01 1.198E-01 2.360E-01 5.066E-01 

1.647 1.246E-01 1.275E-01 2.829E-01 5.351E-01 

1.729 8.349E-02 1.489E-01 3.315E-01 5.639E-01 

1.800 2.790E-02 1.868E-01 3.796E-01 5.943E-01 

1.900 2.390E-02 1.824E-01 4.073E-01 6.136E-01 

2.000 1.870E-02 1.823E-01 4.433E-01 6.443E-01 

2.200 1.330E-02 1.851E-01 4.784E-01 6.768E-01 

2.500 l.OlOE-02 1.939E-01 4.713E-01 6.753E-01 

3.000 8.700E-03 2.202E-01 4.572E-01 6.862E-01 

3.500 9.100E-03 2.398E-01 2.872E-01 5.361E-01 

4.000 l.OOOE-02 2.497E-01 2.324E-01 4.922E-01 

4.500 1.120E-02 2.638E-01 1.937E-01 4.688E-01 

5.000 1.260E-02 2.797E-01 2.389E-01 5.312E-01 



